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ABSTRACT 


A model of a TWS radar is developed that provides a realistic 
computer simulation for comparing various radar tracking methods. 

Prediction accuracy of a simplified alpha - beta tracker is 
compared to that of an adaptive filter. In addition, the effect 


on radar tracking of a variable gate size correlation technique 


is investigated. 
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CHAPTER T 
INTRODUCTION 

A need exists for more extensive use of computer simulations 
as an aid for defining and obtaining the solutions to systems engin- 
eering problems. To gain some insight and experience in this field 
a computer simulation of a Шинээ есені point s a system 
is carried out. 

The criteria for defining the radar parameters and method of 
operation were derived from the following situation. 

A к PM. ship should be able to provide itself with 
_adequate protection against a multiple aircraft raid while transiting 
singly. In order to provide timely detection and adequate нө ese 
for launch of a defensive weapon fue к radar must be capable 
T the following functions: 

Full hemispherical coverage. 
Medium range search i.e. 50-100 n.m. 


Multiple target track while scan. 
High.data rates in critical threat sectors. 


= (00 А Ы 


In Bassin: to minimize costs the radar design should be compa- 
tible with existing weapons. To provide & realistic approach, only 
those designs and techniques capable of being implemented with curren- 
tly available components were allowed. 

A phased array antenna with a monopulse beam was chosen to meet 
the requirement of high data rates on multiple targets. Phased arrays 
are inherently limited in scan. To meet the requirement of full 
360° coverage it was necessary to rotate the antenna mechanically in 
azimuth. Formulation of the problem in this way led to the selection 


of two modes of operation. The first is characterized by complete 


radar coverage with mechanical rotation in azimuth and electronic 
scanning in elevation, with monopulse in elevation angle only. This 
mode has limited tracking capability due to a relatively low revisit 
rate. Mode 2, a sector scan with an electronically steered beam in 
azimuth and elevation and monopulsing in two coordinates, is capable 
of providing 0° to 90° elevation coverage and 3 45° in relative а2- 
imuth. 

It is felt that the two-mode solution is justified since one can 
assume with some confidence that the raid will be radially inbound 
with little or no intentional maneuvering at ranges greater than twenty- 
five nautical miles. Once the raid has been detected and evaluated as 
such by the observer, Mode 2 can be selected to provide a high data 
rate to track the target, whose probability of maneuver increases as 
the range to the target decreases and as his weapon launch point is 
approached. 

The switching range is not critical and can be assumed to be a 
function of the maximum range of the expected weapon available to 
the enemy aircraft (of the order of 5 to 15 n.m.) plus some safety 
time. 

V Once the radar was modeled and the computer program written, the 
Simulation became a convenient vehicle for investigation of the radar 
tracking equations which are a critical facet of a TWS system. A 
comparison of prediction accuracy and resolution of target tracks 
between a classical alpha/beta tracker and an adaptive technique was 
pursued as an extension of the basic premise of utilizing the computer 


as a device for design and analysis. 
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СНАРТЕК ІІ 


RADAR MODEL 


A discussion of some assumptions which affected the choice of 


Mme cype Of radar to be modeled was presented in the introduction. 


This chapter presents some of the reasoning and methods used to select 


basic parameters for the radar necessary to establish a reasonable 


time scale for che Simulation. This initial seetion is followed 


by & discussion of & few of the more importan; subroutines in enough 


detail to supplement the flow diagrams of appendix II where necessary. 


Choice of Basic Parameters. 


Assume: 


(1) 
(2) 
(3) 
(4) 
(5) 
(6) 


ALT (FEET) 


2500 
2000 
1500 


1000 


Target velocity = 600n.m./hr. 

Mittack Yaroiile shown in fig. 

Target weapon release point = 10 m.m. шах. 

Defensive weapon = missile (mach 3). 

Missile requires a high data rate for command guidance. 


Antenna height (h,) = 50 ft. 


TARGET ASSUMED TO BE ON 
“ A RADIALLY INBD. COURSE. 





10 15 20 25 
(Nautical Miles) 


Fig. 1 ATTACK PROFILE 


ШЕЛ) 


Desiring a 30 second sector scan in which to fire and control 
the defensive missle, we must commence sector scan at a minimum range 
of 20 n. m. This ensures an impact point of approximately 15 n. πι., 
providing adequate time for a second attempt before the attacking 
aircraft reaches his weapon launch point. Arbitrarily assuming twice 
this range (MO n. m.) will provide adequate time for threat evaluation. 
Based on the radar horizon given by the following equation, 50 n. m. 
was chosen as the range at which a probability of detection, equal 
to 90% on a single scan, СОО ед sufficient initial detection 


capability. 


r= 21 + \ 2h¢ (miles) 


ht = height of target A/C. 


given : ha SORTU: 


21010): 51488 


r = 41 miles h. 


r 2 5h miles πι - 1000 ft. 


Having obtained an idealized Rmax = 45 n. m., we can calculate 


the pulse repetition frequency. 


R(unamb.) =c = 5 n. m. 
21, 
- 3х108 ш/ 7, 
fo = 5x109 amiee 5.4 X 10 n.m. 
2x 15 n.m. | т. 
f, - 1800 cycles/sec. 
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Assuming a fan shaped antenna beam (1° azimuth by 5.5° in 
elevation) and requiring the equivalent of integrating 4 pulses per 


scan we determine the time for one revolution of the antenna in 


mode l. 
1 è h cycles x 3605 - 0.8 СЕС. 
1800 сусДез/ вес. 1° Beam REV. REV. 
Position 


This figure must be multiplied by the number of 5.5° sections in 
elevation that must be covered. Taking into account the target 
profile and realizing the advantages of short revisit rates a 
compromise value of 54.5° coverage in elevation requiring 8 sec. 


total scan time was’ selected. 
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To simulate the continuous motion of the antenna, azimuth is 
divided into 0.02° increments. The following diagram describes 


the antenna motion. 


DIM? 
19 
m — унь. 

| 
| | 
| | 

| 
| | 
| i 
MEE 

ELEVATION IS | 
| ELECTRONICALLY | 
| SCANNED IN 10 (5.45°) 
INCREMENTS 

| 
| ! 
! | 

! 
| | А 0.0291 
| f 

| | 
| 
| | 

| 5,159 


| аш 


Fig. 2. SIMULATED ANTENNA MOTION — MODE 1 
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Turning now to a discussion of mode 2 or full electronic sector 


scan the following assumptions were made: 


(1) Defensive missile has a 20 millisecond command guidance 
requirement repeated at intervals of 1 second. 
(2) Special purpose beam steering computer capable of: 
(a) Calculating 200 beam positions per second. 
(b) superimposing ; COllimation commands on the beam 
position to obtain variable beam shape. 
(c) Beam shape = 2° x 2° pencil beam. 
(3) 5 ms. illumination time per beam position except 25 ms. 
illumination time for beam positions where targets are predicted 


po be present. 


(4) 10 - target capability, maximum. 


The above assumptions lead to the establishment of the following 
table: 
TABLE 1. Search Beam Positions vs. known targets 


Targets Search Beam Positions Remaining Per 
1 Second Interval. 


0 200 
= 190 
Ц 180 
6 0 
8 160 
10 150 


Ш> 


In a Í 459 AZ. by 60° EL. sector we have 15 x 30 = 1350 (2° x 2°) beam 
positions. For the worst case of 10 targets engaged the revisit 
rate to a "same" search position would be 1350 = approximately 9 

150 


seconds. This scan interval is reduced to 6.75 seconds for the case 


of O targets engaged. 


-45° LIMIT 


JE 







+45° LIMIT 


ELEVATION 


POSITIONS 
9| 190 


ә-------- AZIMUTH POSITIONS ---------» 


Fig. 3. BEAM SCANNING SEQUENCE FOR 2 TARGETS ENGAGED. MODE 2. 


The target detection model for mode I and mode II are essentially 
identical and can be described in the following way. If the radar 
beam passes over a target the target range as given by the target 
generator is compared to the values in Table 2 to determine the 
probability of detection for that range. If Random, the result of 
calling the uniform random number generator RAN1, is less than the 
probability of detection, we assume a detection; otherwise a miss 


or noise return. 
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TABLE 2. Probability of detection vs. range. 


Prob. of Det. Target Range 
` 0.22 60.0 
o2. 55.0 
0.6 54.0 
Ше πο 
0.8 52.0 
сыр 51.0 
0.9 50.0 
0.95 49.0 
0.99 8.0 
0.999 47.0 


Having established a radar target the radar range RADRNG is 
computed by the equation RADRNG = TARGRNG + RANERR where RANERR is 
obtained by calling the GAUSSIAN random number generator and using 


the equation: 


RANERR = DEV * 0.3 


We can see from the above equations that the radar range is 
equal to the true range plus an error term with a standard deviation 
БУ ОШОГО ЫТ ОЗ п.ш. 

" It should be noted at this point that all radar measurements 
are computed in essentially the same manner. That is RADRNG, AZRAD 
and ELRAD are computed as noisy observations of the "true" target 


position, as given by the target generation subroutine, TARGEN. 
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CORRELATION: 

To correlate current radar observations with tracks established 
from predictions on previous observations, subroutine corrass was 
devised. Current radar observations are stored in the dius RADRNG , 
AZRAD, ELRAD, While predicted values from the tracker are stored 
Similiarly. Once during each revolution of the antenna each track 
is compared with each observation; for each observation that a track 
agrees with, a one is placed in the corresponding location of the 
C(I,J) correlation matrix. The columns of the correlation matrix 
correspond to tracks and the rows correspond to observations. For 
& track to agree with an observation certain tolerances were assumed: 


range must be within Í X n.m., azimuth and elevation angles must be 


within = Y degrees. 


X and Y are determined as follows: 


(1) Alpha - Beta filter. 


Range gate = 2.5 n.m. 
Azimuth gate = 2.5 n mn: 
Elevation gate = 1.5 n.m. 


(2) Kalman Filter, 


Range gate 2.0 * SQRTF [P(1,1)] n.m. 


Azimuth gate 1.0 * SqRTF [P(3,3)] п.п. 


2.5 * SQRTF [P(5,5)] n.m. 


Elevation gate 


After this initial correlation several observations may correlate 
with one track or several tracks may correlate with one observation. 


In these situations, a set of four association rules are applied to 
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the correlation matrix, C (I,J), to assign observations to the 


proper tracks. 
These rules implemented by subroutine ASSOC, are as follows: 


(1) A track which correlates with several observations rejects 
any observation held in common with another track, if the common 
observation is the only observation correlating with the second track. 

(2) A track correlating with several observations, some of 
which are not held in common with other tracks, rejects observations 
which are held in common with other tracks. 

(3) When several observations correlate with one track, the 
closest observation is associated with that track. 


(4) When several tracks correlate with one observation, the 


observation is associated with the closest track. 


Rules 3 and 4 are based on the value of range since it has the 
smallest yV? for measurement noise. The above correlation and assoc- 


iation routine is essentially the same as that given by [7]. 
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CHAPTER III 


FILTERS 


Smoothing of raw radar reports and prediction of future obser- 


vations, basic functions of an automatic tracking system, can be 


performed conveniently by a set of equations implemented on a com- 


puter. 


Tn this simulation two distinct schemes were chosen to 


perform these functions. In the simplest case target tracks are 


based on smoothing and prediction of an alpha-beta tracker operating 


in a cartesian coordinate reference frame. 


The smoothing equation is: 


Леп ® Дра t X (An ZA pn? 


The prediction equations: 


Where: 


Ne 


Apa 


An 


igs 


ra~ 


NE =A nel BK, ΓΛ ο.) 


Λίρα “Аза 5 Да %7 


‘Noisy observation of X,Y or Z for the n 


Smooth value of X,Y or Z for the nth scart. 


- Predicted value of X,Y or Z for the Ath scan. 


th scan. 


Smooth prediction of X,Y or Z component of velocity . 


for the nth Scan. 


= Predicted X,Y or Z coordinate for the N + 1 scan. 


T = Time between looks which is essentially constant. 
T = 1 second for sector scan. 
Т = Ө seconds for full scan. 
ои = The smoothing parameter. 
= ae. Сот орала filtering. 
2-0 


Normally associated with the alpha-beta smoothing prediction 
equations are a set of rules for determining bin size. The rules 
are usually based on track firmness, length of time since last 
correlation and radar range. Basically this set of rules for 
varying the correlation gate size is an attempt to make the system 
adaptive ee dynamics. 

Experimentation with the simple alpha-beta tracker indicated 
that given a reasonable fixed gate size based only on accuracy of the 
radar measurements, tracks could be maintained on targets with mild 
maneuvers, ie., less than 3°/sec. turns, with little difficulty given 
a suitable value of alpha. An important result observed from these 
simulations was that the real problem associated with tracking is 
the time delay before sensing a target's maneuver initiation. 

Investigation of a criterion for determining target maneuvers 
seemed to be the next logical step. Changes in the magnitude and 
angle of the target velocity vector in the X - Y plane were looked 
at briefly as one possibility. From the data accumulated, inciden- 
tal to other tasks being performed, the changes in angle appeared 
to be an unlikely candidate due to extreme noise caused by measure- 
ment errors of the radar. The magnitude of the vector, on the other 


hand, appeared to be rather insensitive to target maneuvers. 
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In lieu of pursuing this somewhat uncertain investigation, it was 
decided to use techniques that were in comparison well defined. The 
Kalman filter was chosen as the best solution to the problem for the 


following reasons: 


(1) While the Kalman filter can not predict target maneuvers, 
the P matrix generated by the filter equations is directly influenced 
by target maneuverability (Q matrix) and at the same time is a 
measure of our confidence in predicted values. It seemed reasonable, 
therefore, to make the correlation gate size proportional to the 
corresponding values of the P matrix. 

(2) If an observation is missed due to a target maneuver we 
can expand the gate size to some larger value to increase the pro- 
bability of correlation on the next scan with automatic gate size 


reduction if we succeed. 


The Kalman Filter Equations as given in (4): 


5 t Ë -1 

Ge - PEU {ΠΡ 027 + R] 
А, ^ A 
п 7 "n/n-i * 8n [4 - HXn/n-1] 
Pa/n = Ph/n-1 - Sn#Py /n-1 
m 

ntl/n ` 0 п/п 

Р = ØP 0 +9 

n+1/n n/n 


22 


These equations represented in block diagram form appear in 


mim cd. 


DELAY 


m 
62 
52 
S 


n ; 


Xn/n- 


Fig. 4. BLOCK DIAGRAM OF FILTER EQUATIONS. 
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where: L= HK + V 


Radar Range 
Radar Azimuth 
Radar Elevation. 


N 
li 


or 


lj 


V^ represents the noise associated with the radar measurement 


of the three observables and is gaussian. 


H = observability matrix. 
πο ο ο ο 0 
or Н-1001000 
шт о 


@ = The state transition matrix that represents the dynamics 


of the aircraft target. Assuming the aircraft can be rep- 


i 
resented as pure inertia of a -о plant. 
5 


2з 


Ме have 





Yields: 


hence: 


where: 


we know: 


hence: 


' Assuming the 


0(%) 


x, (0) 


X, (0) 


с(5) 


ЖЕ) х.(8) Хү(8) Хү (8) 


Fig. 5. FLOW GRAPH OF 1/S^ PLANT. 





X (8) Su (S) + 1 Χο (8) + 0 U(S) 
X (s) Ξ 0 x (8) 40 e (5) + 1 U(S) 
fo 1' INS 
A = => Ότι... 
αμα ο 1/5 
1 
Δ = σα 
5 1 3 р : 1/8 1/54 
[opem zL159 πι. = |0 1/8 
А А 
1/8 1/82 
ф (тэ ота 010 τος 
1 t 
0,(6) 5 
O 1 


components of the X vector are uncoupled. 


SRO OLOO 
© © cr 
ссу оте ого 
ΓΞ 


on © оғо © 
сас CO OOo Oo 


t 


1 


The covariance of the measurement noise is by defini- 


t 


tion = F [V V ] 


In this case the- standard deviation of the measurement 
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noise on Range, Azimuth and elevation are Known to have 


the following values. 


Range = 0.3n.m. 
Az = 0.5 deg 
El = 0.5 deg. 
hence: 0.09 0 0 
R = ог. 0 
0 Ш 0225 


X n/n-l represents the predicted states of & radar target. 
To initialgswewbhewRilter for a^particular target, that is specify 
X 1/0, the following method is used. 


Radar Range 


к = шо 
Х 1/0 = Radar Azimuth 
AV - 0.0 


Radar Elevation 
ВІ 0.0 


This method assumes the target is closing the radar radially 
(А7 - 0.0) аба Қ altitude (El © 0.0) with a velocity of approx- 
imately 600 n.m. /hr. 

Р п/п-1 represents the amount of uncertainty in the predicted 


values of target states X n/n - 1. 


where: is the true state. 


Х 
^ 
Және іле песе value. 


To initialize the filter we must provide. 


Eoas Given" xX defined previously and R the covariance of the meas- 


1/0 10 
urement noise; 


Р (1,1), Р (3,3), Р (5,5) аге Кпомп. 


Again relying on the assumption of a radially closing target at low 


go 


altitude with a speed of 600 m.m./hr. we can assume the remain- 
2 2 e 
ing diagonal elements of P which represent Ve: Vie , NE 
will be small. Actual values were determined empirically and as а 


result the initial covariance of error matrix is determined. 


0.09 0 0 0 0 0 
0 0.0278 0 0 0 0 
Р- 0 0 го 0 0 0 
0 0 0 ΠΕ» 0 0 
0 0 0 0 0725 O 
O 0 0 0 0 0.00109 


Th нда be noted here that Ῥ] /0{ 333) = 2.5 is used in place 
of 0.25 to aid tracking maneuvering targets. If an observation fails 
to correlate with an established track the P matrix for that track 
is set equal to the above values. The gate sizes, and in particular, 
the Azimuth gate, are directly related to the P matrix values. This 
procedure enlarges the gate size thus 30275771: the probability of 
correlation on the next scan. 

Lastly but certainly not of the least importance is the Q matrix 
specification. The Q matrix which represents perturbations to the 
state vector due to target maneuvers must be representative of the 
full range of movements the target is expected to use. 

Values for Q were “ог by -——Ó "worst" case such as 
6000 ft. /min. rate of climb, 1 g. linear acceleration and 12° / 
second turning rates. One half these "worst" case values were 
actually used in an attempt to average the Q matrix over the full 


range of expected maneuvers. 
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Estimation of Q (1,1) 
Assume the aircraft can vary speed by (100 n.m / hr. max.) 
ο ο Dc n m. /see. 


А 


Q 


Е [w w°] where W is a random 

Signal caused by target maneuvers. 

Taking 1/2 (0.027) implies 

Eg» = UNUS ο ο - 0.0169 

Estimation of Q(2,2) 
Assume the aircraft can accelerate at (1g) max. 
тыйы Ор" (їй = .o027 nm., х 1 sec. 


See 
imprres thay @la.e) = (.0027) (.0027) 


гм 


. 00001 
Estimation of Q(3,3) 
Assume the aircraft can make & (12? / second) max. turn. 


Taking 1/2 (12°/sec) = 6°/sec. and using fig. 6 


We have: 
‘chord length = 2r sin 9/2 
21:26 зө 
30. 
Непсе. 


0 (3,3) - (.083) 4088) - .00692 
Estimates for the remaining terms of Q were obtained in a like 


manner with Q taking the following form. 


0.0169 0 0 О 0 0 
0 200001 0 0 0 0 
0.0060 0 0 0 
Q est.'z 
0 0 О ο О О 
0 0 0 0 ο ο] 0 
0 0 0 0 0 ο ο! 


СТ 


EO 
X = RAD Oc we On 
qoe 12082229 
Vel = 600 m ma /h = 0 o 00010100 
η: 
: „Л 
| jag 
DE 
: 
i tN 5 
è] e 
= 
| Е 
Hé- 
| 


"ig. Ὁ. CHANGE IN RADAR AZIMUTH DUE TO AÀ TURNING TARGET. 


Simulations were run to optimize the Q values in order to obtain 
the best filter response over the range of 0°, 3°, 6°, and 12° /sec. 


turns with constant speed targets. 


0.00169 0 ο 0 0 ОЕ 
0 0.0030 0 0 0 0 
0 ο. 0.015 0 0 0 
Q exp. = | 0 0 0 Ὁ. 0 0 
0 0 0 0 0.0001 0 
0 0 0 0 0 0.00001 


Тһе Qe matrix values where not increased in magnitude as might 


at first be expected but were set equal to Qi for nhe RolNowinP reasons: 


(1) The eight second scan time applies to targets of fairly 
long range compared to target ranges anticipated for sector scan, and 
as such the targets are expected to maneuver as little as possible 


in an attempt to close the ship. 


(2) It is intuitively felt that T(ave.) for aircraft maneuvers 
for this type of target might be expected to be 3-5 seconds vice 8 


seconds. 


(3) The effect of maneuvers acting on a spherical coordinate 


system decreases with increasing range. 
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CHAPTER IV 
CONCLUSIONS 

The radar system simulation presented in this paper is felt over- 
all to be & simplified but realistic model of & sophisticated system 
which could be implemented with current "state of the art" hardware. 

The basic main program / subroutine organization divides the 
system into easily recognized functions capable of being extended in 
scope, modified in part, or replaced entirely with & new concept, 
while ιο ρα the integrity of the remainder of бие ресоспай, 

Any portion οἱ ΕΠ; can be reworked to meet special requirememm 
ОТ & particular problem of interest. 

By the addition of new subroutines which might include phenomena 
associated with clutter, target scintilations, and atmospheric atten- 
uation a more realistic detection probability, given a set of in- 
itiad Conditions , сопа, асоче геа 

Hence the original objective of building a model for simulation 
of a hypothetical radar system as a tool for analysis of a variety 
of radar system problems has been met. 

As an exercise with practical significance, the simulation was 
used to investigate the relative effectiveness of several types of 
radar tracking filters. 

-wa While no pretension of an extensive analysis is implied, the 
author feels that a reasonable and unequivocal comparison of the 
filters can be made from the material presented. 

Based on the ensemble averages of the squared difference between 
predicted and actual positions the Kalman filter obviously provides 


a somewhat better tracking response for all target tracks tested. 
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The tracking ability of both filters appears to Ре about equal 
for "look alike" targets in close proximity. The ability of the 
tracking routine to resolve these targets depends on using relevant 
gate sizes. The gate size in turn is directly related to the accu- 
racy of the radar measurement of range, and angles, and to the extent 
of maneuver capability of the target. The minimum gate size allowable 
without inducing excessive non-correlation due to observation noise 
would thus be essentially the same for both filters. 

The most significant — of Kalman Filter tracking appears 
to be in its ability to track maneuvering targets without dropping 
track. This advantage is a result of its ability to automatically 
increase the gate size if & firm track fails to correlate, thus im- 
proving the probability of correlation on the next scan. 

The requirement for peripheral tracking functions such as track 
firmness and quality would surely be reduced in a system using Kalman 
filtering since the (P) matrix can be used as a measure of this type 
of imformation. This feature affects somewhat the disadvantage of 
increased computation time due to the recursive algebraic equations 
required to calculate the G and P matrices. 

Finally, the dependence of the Kalman filter on the Q matrix 
dictates that an efficient and reliable method for estimating Q for 


a particular target is essential. 
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APPENDIX I 


GRAPHICAL OUTPUT 


Graph Page 
Alpha-beta filter 09/sec turn rate. 36 
Kalman filter 09/sec turn rate. 3T 
Alpha-beta filter 39/sec turn rate. 38 
Kalman filter 39/sec turn rate. 39 
Alpha-beta filter 69/sec turn rate. 40 
Kalman filter 6°/sec turn rate. hl 
Alpha-beta filter ТРА rate. 12 
Kalman filter 12°/sec turn rate. 43 
Alpha-beta filter Crossing targets T= 8 «08. hh 
Kalman filter Crossing targets T = 8 sec. 15 
Alpha-beta filter Maneuvering target T = 8 вес. 16 
Kalman filter Maneuvering target T= 8 sec. 7 
Alpha-beta filter Crossing targets T= 1 sec. 48 
Kalman filter Crossing targets Т = 1 sec. 49 
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Interpretation and constraints 
(1) The first eight graphs provide a comparison of the pre- 


diction accuracy of each of the two filters investigated. 


The following rules apply: 
(a) Four different target tracks are employed i.e. 0°, 

35 “6° and 12°77 seer | 

(b) Target turning rates are initiated at T = 6 sec. 

(с) Each graph represents the results in range (m), 
azimuth (X), and elevation (4) of a "Monte Carlo" ensemble average 
oí ПУСК Ла, 

(4) Sample rate = 1 sec. appropriate to Mode 2 operation. 

' 


(e) Alpha-beta filtering is performed in a cartesian coordinate 


reference frame while Kalman filtering is done in spherical coordinates. 


Filtering of the azimuth coordinate, which happens to be changing 
most rapidly, is about equal for both filters. The Kalman filter 
is obviously superior in range and elevation predictions. 
The following 2) table applies for all succeeding graphs. 
true track 
Х predicted position target one. 
А predicted position target two. 
(2) Graphs nine and ten represent filtering ability for "look 
alike" targets. i.e. Же range, elevation, and azimuth. 
The ability of both filters, supplemented by a set of correlation 
and association rules, to distinguish between similar targets in close 


proximity leaves something to be desired. A reduced revisit rate i.e. 


(1 sec vice 8 sec.) has little or no effect. (see graphs 13 and 14). 


3h 


It seems reasonable therefor to assume that a higher resolution radar 
is the only solution to this, problem. 

(3) Graphs ten and eleven indicate the improved tracking response 
of the Kalman filter for maneuvering targets, aided by & variable gate 
size technique. Investigation of the printed output for this run shows 
that the Kalman filter operating on spherical coordinates consistantly 
predicted the position of the turning target within the limits of the 
gate size. | 

In contrast the Alpha-beta filter was unable to follow azimuth 


changes with sufficient accuracy to obtain the necessary correlations 


and as a result dropped track excessively. 
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APPENDIX 11 


LIST OF VARIABLES 


This appendix contains a listing of the variables used in the 


simulations, arranged in alphabetical order. 


Variable 
ALPHA 
ALT(K) 
ALT1 
AZBEAML 
AZBEAMT 
AZDIFF 
AZDOT 


AZL (K) 


AZG(K) 
AZ2(T) 
AZRAD(K) 
μπα 
PK) 
AY 
BETA 


BOXMAX (I) 
BOXMIN(I) 


ο τα) 


СОМТ(К) 


DELTA 


Definition 


Position smoothing constant of the X-B tracker. 


Altitude of the K"? target. 


A dummy variable. 

The leading edge of the radar beam. 
The trailing edge of the radar beam. 
Predicted azimuth - true azimuth. 
Azimuth rate. 


Counter for the number of times the «th 


fails to correlate. 1 


Azimuth gate size of the гн 


track 


track. 
The ensemble average of ELDIFF %%2 
The noisy radar azimuth of the geo target. 


A dummy variable. 


- ™% 


The true azimuth of' the ΚΟΠ 


target. 
A dummy variable. 


Smoothing constant for velocity in thea-B tracker 


The maximum value allowed for a radar return if 
it is to be associated with a particular range box. 


The minimum -value allowed for a radar return if 
it is to be. associated with a particular range box. 


Observation / track correlation matrix. 


Azimuth of tHe leading edge of the radar beam on 
initial detection of the Kt*h target. 


Normalized velocity vector. 
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Variable 


DOTT 


DOTX 


DOTY 
DOTZ 
ЕГ,ВЕАМВ 
ELBEAMT 
ELDIFF 
ELDOT 


ELERR 


ELG(K) 
ELRAD(K) 
ELT(K) 


ELTEST 


FLAG 
G(I,J) 
GO 

HOT J) 
HD(K) 
HPHT 
HT(I,J) 
IA 
IAT(I) 


IAZ(I) 


IAZREF 


Definition 
A dummy variable. 


The predicted velocity components of the Ke target 


as given by the(-B tracker. 


Bottom extemity of the radar beam. 
Top extemity of the radar beam. 
Predicted elevation - E elevation. 
Elevation angle rate. 


Radar measurement noise added to the true elevation 
angle. : 


Elevation gate size of the qth track. 
Noisy radar value of target elevation.’ 
True target elevation angle. 


Parameter for distinguishing between range boxes by 
elevation angle. 


Flag is set by MOVAVE if a target is detected. 
Gain matrix for the Kalman filter equations. 
Used by SCAN 360 to compress full scan. 
Observability matrix for Kalman filter equations. 
Heading of the коп aircraft in degrees. 

A dummy variable. 

The transpose of the observability matrix. 

Range box index. 

Used as a counter in conjunction with IAZ. 

Radar return signal counter which is compared to 
a threshold value to determine whether or not a 


detection has occured. 


kefercnce azimuth of sector scan. 


JL 


Variable 
TI 
INDE(I) 


INDEX (I) 


ITEL 
ITELI 
ІТЕМР 
ІТЕМРІ 


KO(I) 


KER 
LIMB 
LIML 
LIMR. 
LIMT 
NB 
NBT 


NO 


NR(I) 


ΝΤ 


NUM1 
NUM2 
NUM3 
NUMSCAN : 


NUMTRK 


Definition 


Elevation beam position indicator. 

SEE INDEX (I). 

INDEX in conjunction with INDE determine when a track 
fails, two consecutive times, to correlate thus in- 
itiating drop track prodedure. 

The fixed point form of TEL. 

A dummy variable. 

The fixed point form of TEMP. 


A dummy variable. 


A flag used to initialize the Kalman filter for a 
particular observation. 


Flag to indicate an improper matrix'inverse operation. 
Lower limit of the sector scan. 

Port limit of the sector scan. 

Starboard limit of the sector scan. 

Upper limit of the sector scan. 

The number of runs. 

Number of discrete times. 


Counter for the number .of target bins currently 
active. 


Flag to indicate known targets in sector scan. 
Number of targets. 


Sector scan counter. 


Dummy variables. 


The number of 360 degree scans of the radar. 


Number of current active tracks. 


De 


Variable 


OISE 


P(I,J) 
PAZRAD 
PAZY 


PD(I) 


PDET 
PELRAD 


ΕΠΙ τα”) 


BEINNF(T, J ) 
PHI8(I,J) 
FHIT(T,J) 


PTTCH(K) 


PITCH1 
ποτ) 
PRADRNG 


PRED 


PXRAD 
PYRAD 
PZRAD 
BUT. J) 


Q1(I,J) 


8 (1.2) 


R(I,J) 


Definition 


A dummy variable for the noisy radar observations 
(Range, Azimuth or Elevation. ) 


The covariance of error matrix. 
The predicted radar azimuth of the target. 
A dummy variable. 


A discrete value of probability from the proba- 
Баша ужо κος ος Table. 


Præsability of deection. 
Predicted value of the target elevation. 


The state transition matrix used in the Kalman 
filter. 


The state transition matrix for the one second case. 
The state transition matrix for the eight second case. 
Transpose of the PHI matrix. 


Value of the nose attitude of the кеп 


сонс 


target in 


A dummy variable. 

Storage for P matrix associated with each track. 
Predicted radar range of a target. 

A dummy variable used in the Alpha-Beta tracker 
for the predicted value of range, azimuth or ele- 


vation. 


Predicted X,Y,72, components of target position. 


Matrix representing target maneuver capability. 


Q may be set equal to QI or Q8 as required by the 
appropriate scan time. 


Covariance matrix of the observation noise. 
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Variable 
RADDEG 
RADIFF 
RANDOM 
RANERR 


RANGE(T) 


RADRNG(K) 
RADRNGL 
RDOT 


RG(K) 


S(I) 


SMOOTH 


SUM 


ШІ 
T2 


TARGRNG 


TEL 
TEST 
TESTT 
TNT 


THETA 


VEL(K) 


VMAG2 (I) 


Definition 

Constant used to convert radians to degrees. 
Predicted range - true range. 

n distributed random number. 


Error in measuring target range. 


Discrete value of range used to enter the pro- 
bability of detection table. 


Noisy target range as given by the radar. 
A dummy variable. 
Range rate. 


Range gate size for the m track of the cor- 
relation routine. 


Α dummy variable used for temporary storage. 


The filtered radar observation in range, azimuth 
or elevation. 


Temporary value of the radar beam leading edge. 


Time interval for updating targets equal to scan 
time. 


Dummy variables. 


True value of target range as given by the target 
generator. 


Temporary sector scan beam reference. 


Value of target heading in radians with refer- 
ence to the positive X-axis. 


Speed of the Ken target in n.m./hr. 


Ensemble average of radiff. 
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Variable 


VNB 
VPHI2(I) 
ШЕ Т.) 
ХРІ 


XRAD(K) 


XRAD1(K) 
XRAD2(K) 
XS 

XT(K) 
XT1(K) 
XT2(K) 
YRAD(K) 
YRAD1(K) 
YRAD2(K) 
ИРК) 
ҮЛЕК) 
YT2(K) 


207). 


ZRAD(K) 


ATEK) 


Definition 


Fixed point form of NB. 
Ensemble average of azdiff. 
Predicted state vector. 

A dummy variable. 


X coordinate of the κο 


radar. 


target as given by the 


Stored radar track quantity for graphical output. 


Smooth state vector. 
Шы; : th 
The true X coordinate of the K target 


Stored target track quantity for graphical output. 


Radar Y coordinate of the oe target. 


Stored radar track quantity for graphical output. 


тетен жесш тате Өте κ target. 


Stored target track quantity for graphical output. 


Measurable quantities input to the Kalman filter 
equations. 


The radar Z coordinate of the КОЗ target. 


The true Z coordinate of the Kee target. 


1; 


Note: Variables associated with the following input/output апа 
‘matrix algebra subroutines, have not been defined, since 
they are essentially dummy variables. 


1. ΑΠΡ 
2. PRINTT 
3. PROD 
|, READD 
5. КЕСІР 
6. TRANS 
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This Appendix contains a 


below. 


APPENDIX III 


LOGIC FLOW DIAGRAMS 


series of Logic Flow diagrams as listed 


MAIN PROGRAM 


SUBROUTINE 


SUBROUTINE 


SUBROUTINE 


SUBROUTINE 


SUBROUTINE 


SUBROUTINE 


SUBROUTINE 


SUBROUTINE 


SUBROUTINE 


SUBROUTINE 


SUBROUTINE 


2T 


TARGEN 
SCAN 360 
DETECT 
MOVAVE 
CORRASS 
ASSOC 
ABFILT 
KALFILT 
SETUP 
TWS 


MONTE 


LIO CIO SL 


FLOW CHART SYMBOLS 


A connector or terminal. 


An offpage connector. 


A predefined process or module/subroutine. 
A more detailed flow chart of this subroutine 
is also included. 


Input/output other than display. 


Decision. 


Processing, annotation. 
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MA\N PROGRAM 







START 


CONSTANTS 
FLAGS 









READ * of TARG. 
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END 
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VARGEN 


START 


«5 
D——| HD: HU * AuD 


NO V 


Yes 
a VELZVEL+AVEL 
| 
| 


XES 
T PITCH: Pe ΝΟ... 
U 


OP 


NO 





ITEXT p o 
YT aka +: 
TaT e 













RECTANGULAR 
TO SPHERICAL 
CONVERSION 





ALL 
TARGETS 
UPDATED 
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SCAN 360 
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AZIMUT IA < ooo 
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ELEVATION 





MODIFY 
AS 
No REQUIRED 
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+ 0.3% DEV 
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65 







<> 
- NO 
ν 


ti 
М 
y 
RETURN 
V 


No 








Az4 (obs) 
= А? Ҷ(0ь5у 
-1.0 





REJECT 
COMMON OBS. 


V 


NO 


ONT) 
ου (С 
ASS 





Е 
IAT 
soc 

с озеет 


PAIR 







SSOCIATE 
А 


т 
CLOSES 





67 


ABFILT 


(50:1с-0 
TRACK NO. =L 
СЭР” 


Y 
δος. NO. = ҮС 
U 
° » ҮЕ5 


SPHERICAL 
To RECT. 
CONVERSION 
























UPDATE THIS 
TRACK BY 
ASSUMING 
OBS. = PRED CL) 












SMOOTH 
PRE Dict 
x x 


SMOOTH 4 
PRE DICT 
= 





KALFILT 


START 





READ HRÓQ 
CALL. TRANS 


T φῖ 














CALCULATE 
BEQ Gies ΙΝ 
EILTER CHAP. 
CALL PROD, ADD, 
TRANS, RECIP 


CALCULATE 


G 


CALCULATE 


Xs 












CALCULATE 


Puin 


CALCULATE 
Xu Ig -1 


CALCULATE 
P % [na 












INITIALIZE 
KC (y=1---10) =O 
РКІ т Рамат 








YP YZ RA0RNGGO 








X PI = PRADRNG (L) ХРА = PRADRACO) 


1 
i 






; : 
Р -. РКОО 


Р = Ῥινιτιαι. 
Ἔ = КАЎКА (το 


Z= PRK0RMG (D 


$ 
P< Рутул 
2>казуама(се) 
L = 







581 Ц р (сомт) 


Νο 


RADRNG(LS Ξχδι 


b 
i 


PRADRNG(L) = Хе! 


, 
t 
% 


RG(L) = K «Sart 
(Р (1,05) 









сек 












INITIALIZE 
PRINT OUT ARRAY 
K= RAID TARG. NO. 
JAZ REF = AZRADIK 
NRCK) = 











SET LIMITS OF 
SECTOR SCAN- 
МАХ ВЕ МоО1- 
FIED PRIOR TO 
RUN TIME 


V 
WUMB -мим8-4-1 
/ 





TRACKS 
CWECKED 


DROP 
TRACK 












SET BEAM 
POSITION 
FOR KNOWN 


TARGE T 


ЇЭ 


p 
4 
C 


NO 
RADRNG > O 


ος 


< 


ES 


AZRAD =AZT+4 ERR 


ELRAD=ELT + ERR 


1! 


ғ) 


τ 


с 


FLAG3=0.0 


ING = 3% of BEAM 
POSITIONS 


SET [th 
BEAM 
POSITION 


NXT 
TARGET 





YES 


i 


TWS 


V 


APRAD =AZT+ERR C 


М 


ELRAD * ELT *ERR 


ү 
U 














BEAM 
POSITION 
= INC 


NO 


YES 





CALL 
TARGE N 
ASSOC 
Альт 


CALL ү 
CORRASS 
SETUP OUTPUT 
u 


V 


€ 


T По” FLAGS = FLAGS «1.6 
| An 


NUMTRK= NUM-Ke| 
NR (<) = K 










e 





76 
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